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We propose a new method for evaluating incomplete and complete fusion cross sections 
separately using the Continuum-Discretized Coupled-Channels method. This method is 
applied to analysis of the deuteron induced reaction on a 7 Li target up to 50 MeV of the 
deuteron incident energy. Effects of deuteron breakup on this reaction are explicitly taken 
into account. Results of the method are compared with those of the Glauber model, and 
the difference between the two is discussed. It is found that the energy dependence of the 
incomplete fusion cross sections obtained by the present calculation is almost the same as 
that obtained by the Glauber model, while for the complete fusion cross section, the two 
models give markedly different energy dependence. We show also that a prescription for 
evaluating incomplete fusion cross sections proposed in a previous study gives much smaller 
result than an experimental value. 

The understanding of the fusion reaction mechanism is one of the most impor- 
tant and challenging subjects of nuclear physics. Description of incomplete fusion 
processes, in which a part of the projectile is absorbed by the target nucleus, with 
emitting other projectile fragment(s), is particularly interesting and important. So 
far, some theoretical models of the incomplete fusion, also called breakup fusion or 
inclusive breakup, of a two-body projectile have been proposed. In these models, 
the incomplete fusion reaction was described as two-step processes, i.e., the projectile 
is broken up first and then one of the two constituents is absorbed by the target. The 
calculations of the fusion cross sections were carried out by using the Distorted Wave 
Born Approximation (DWBA) assuming that the emitted fragment can be treated 
as a spectator in the final state. Recently, roles of breakup (continuum) states of 
a weakly-bound projectile in the fusion reaction have been discussed 4 )' 5 ) using the 
Continuum-Discretized Coupled-Channels method (CDCC). 6 ) CDCC was proposed 
and developed by Kyushu group, and has been successfully applied to analyze vari- 
ous reaction processes; see, e.g., Refs. 7)— 9). In Ref. 4), an attempt to calculate the 
incomplete and complete fusion cross sections separately with CDCC was described. 
The assumption used in the separation of the two was, however, unrealistic for some 
reasons; we will return to this point later. 

The incomplete fusion process in a deuteron induced reaction on Li targets at 
incident energies up to 50 MeV attracts wide interests of not only nuclear physi- 
cists but also nuclear engineers, because the emitted neutrons through this process 



E-mail: hashimoto.shintaro@jaea.go.jp 



2 



Letters 



are planning to be used in the International Fusion Materials Irradiation Facility 
(IFMIF). 10 ) Understanding of the reaction mechanism of this incomplete fusion pro- 
cess, or, equivalently, the inclusive (d, n) process, with evaluating theoretically the 
absolute value of the cross section is necessary. Moreover, nuclear data of inclusive 
(d, n) reactions on other various targets such as Be, Ta, and W are of crucial im- 
portance for studies on accelerator-based applications, i.e., shielding of the deuteron 
accelerators including IFMIF, and medical applications for Boron Neutron Capture 
Therapy (BNCT). Very recently, Ye et al. 11 ^ showed that the main part of the dou- 
ble differential cross section (DDX) data 12 ' of the emitted neutron from the inclusive 
7 Li(d, n) reaction at 40 MeV is reproduced very well by the proton stripping cross 
section d 2 o^r/ \dE n df2 n ) , which corresponds to the (d,n) incomplete fusion process 
in our terminology, added by the elastic breakup cross section d 2 aEB/(dE n df2 n ); E n 
and fi n are the energy and solid angle of the outgoing neutron. In their study, 
d 2 a^ R /(dE n df2 n ) and d 2 cteb / {dE n df} n ) are obtained by the Glauber model and 
CDCC, respectively, and the reason for the surprising success of the Glauber model 
at such low energies (~ 40 MeV) was discussed. 11 ) Nevertheless, it is important 
to evaluate the accuracy of the Glauber model calculation of d 2 a^- R /(dE n df2 n ) be- 
low 40 MeV, in which experimental data are very scarce, by comparing it with that 
obtained fully quantum mechanically. 

In this Letter, we propose a new approach for calculating the complete and in- 
complete fusion cross sections separately using CDCC. As an important advantage 
of the present method to the preceding studies, 4 )' 13 ) we separate the two fusion pro- 
cesses by the physics condition on the absorption of each constituent of the projectile 
by the target nucleus. In our model, a possible contribution of the breakup channels 
to the complete fusion process, as well as that of the elastic channel to the incomplete 
fusion process is properly taken into account. As we mention below, the new method 
contains a free parameter, i.e., the absorption radius. This parameter is determined 
using the result of the Glauber calculation at 40 MeV that can be interpreted as an 
experimental value of the proton-stripping incomplete fusion cross section. We then 
apply this method to the Li(d, n) reactions at different deuteron incident energies 
from 10 MeV to 50 MeV. The (d,n) and (d,p) incomplete fusion cross sections and 
the complete fusion cross section thus evaluated are compared with the results 11 )' 14 ^ 
obtained by the Glauber model. Note that we focus on the cross sections integrated 
over emission energies and angles in this work. 

We describe the 7 Li(<i, n) reaction with the three-body system shown in Fig. [TJ 
R is the relative coordinate between the Li target and the center of mass of d, and 
r is the relative coordinate between p and n. The coordinate of p (n) relative to 7 Li 
is denoted by r p (r n ). The three-body Hamiltonian is given by 

H = T R + U p (r p ) + U n {r n ) + V Conl (R) +T r + V^(r), (1) 

where Tr and T r represent the kinetic energy operators associated with R and r, 
respectively, U p (U n ) is the optical potential between 7 Li and p (n), V Coul is the 
Coulomb interaction between d and 7 Li, and V pn is the interaction between p and 
n. Note that we neglect the Coulomb breakup processes in this study, since we are 



Letters 



3 




Fig. 1. Illustration of the three-body model of p, n and 7 Li. 



interested in the (d, n) reactions enough above the Coulomb barrier energy. 

In CDCC, the three-body wave function \P(R,r) is expanded in terms of the 
eigenfunctions of the p-n system {<Pj(r), i = 0-i max }: 



'max 



JM i=0 



JM 



(2) 



where <Pq represents the ground state of d and <Pi (i ^ 0) the ith discretized continuum 
state. The expansion coefficient Xi{R) describes the relative motion between d in the 
zth state and Li. The discretized continuum state 4>i is obtained by the so-called 
average method 6 ^ as 



#i(r) 



1 



k L 



dk<P(r, k), 



(3) 



where <P(r, k) is the p-n scattering wave function with the relative wave number k, 
and Aki = ki — fej-i- <P(r,k) satisfies 



[T r + V m (r)]$(r,k)=e$(r,k), 



(4) 



where e = h 2 k 2 /(2/j, r ) with fi r being the reduced mass of p and n. 

The three-body Schrodinger equation using the wave function of Eq. (J2J) is given 

by 



(H - E)\P(R, r) =0, 



(5) 



where E is the total energy. Multiplying Eq. ([SJ by <P* from the left, and integrating 
over r, we obtain the following coupled-channel equations for Xi(R) : 



(t r + V p c ™\R) +e l -E) Xi (R) =-J2 FdR)Xj(R), 



(6) 



where £j is the internal energy of the p-n system in the ith state and 

F^R) = (SjKUp + Un^r 

is the coupling form factor. Equations (|6|) are solved under the usual boundary 
conditions for Xi(R)- 6 ^ 
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The imaginary part of the optical potential is considered to describe the particle 
absorption by the target nucleus. Thus, the fusion cross section (absorption cross 
section) is given as the expectation value of the imaginary part with the wave function 
of the system. This cross section contains both contributions from the complete and 
incomplete fusion processes; we henceforth call this the total fusion cross section 
o"TF- In the present three-body model calculation, cjtf is obtained by 

VTF = ^\(*\(-W P -W n )\V)\, (7) 

where W p iW n ) is the imaginary part of U p (U n ), hr is the reduced mass between d 
and 7 Li, and Kq is the d- 7 Li relative wave number in the incident channel. Note that 
the integrand on the right hand side (r.h.s.) of Eq. ([7]) is compact (L 2 integrable), 
since we discretize the p-n scattering wave functions with Eq. ([3]). Another important 
point to be noted is that the imaginary part of the nucleon- Li optical potential 
describes not only nucleon absorption but also other processes such as the inelastic 
scattering to the excited states of Li. Since the nucleon absorption has the main 
contribution to the r.h.s. of Eq. ([7J), however, we regard it as the total "fusion" cross 
section as in many other studies on fusion reactions. 

To separate the (d,p) and (d, n) incomplete fusion cross sections, a[p and crjp^ 
respectively, from otFi we divide the integration region in Eq. (|7|) as follows. The 
explicit form of the expectation value on the r.h.s. of Eq. ([7]) is given by 

\{#\(-W p -W n )\9)\=-J dr p J dr n V*(R,r){W p (r p ) + W n (r n )}V(R,r), (8) 

where we have changed the integration variables from (R, r) to (r p , r n ). We now 
separate the integration regions over r p and r n as 

dr p / dr n = dr p dr n + dr p dr n 

J Jr p <r* h Jr n <r^ Jr p <r^ h Jr„>rJ b 

+ / dr p / dr n + / dr p / dr n , (9) 



•/r p >r£ b Jr„<r% h J r p >r p lh Jr n >r^° 

where r^ b (c = p or n) is the interaction range of W c ; at r c > r^ b , W c is assumed to 
be negligible. The first term on the r.h.s. of Eq. © corresponds to the process in 
which both p and n are located within the range of W c and absorbed by 7 Li. In the 
second term, p is assumed to be within the range of the absorbing potential, while n 
is free of the absorption. Thus, it gives the integration region corresponding to the 
(d, n) incomplete fusion process. Similarly, the third term corresponds to the (d,p) 
incomplete fusion process. It is obvious from the definition of r^ b that the fourth 
term gives no contribution to cttf- Schematic illustration of these four integration 
regions is shown in Fig. [2j Using Eq. Q , otf is decomposed into the complete fusion 
cross section <tcf and the above-mentioned two incomplete fusion cross sections, i.e., 

are = *CF + <T$ + , (10) 
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Fig. 2. Schematic illustration of the four integration regions. See text for details. 



where 

CCF 



dr r 



r v <ri, 



dr n V*(r p ,r n ){W p (r p ) + W n (r n )}V(r p ,r n ), (11) 
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r„<ri 



dr n &*(r p , r n )W p (r p )V(r p , r n ), (12) 



dr n V*(r p ,r n )W n (r n )&(r p ,r n ). (13) 



The total wave function ^(r p ,r n ) is obtained from &(R, r), which is given by the 
CDCC calculation, by the straightforward transformation of the variables. 

We remark here that the above expressions of the three components of the total 
fusion cross section are obtained by properly considering the physics condition on 
the absorption corresponding to each process as mentioned above. On the other 
hand, following the definition of the incomplete fusion cross sections of Refs. 4) and 
13), CTjp (c = p or n) could be expressed by 
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(«0 



r) 'max 

~ Z ^R \ -» \ -> 
IF.prev ~~ fc2 

11 ^° JM i^Q 



[urn 



JM 



W r 



JM 



(14) 



i.e., the integration was done in the entire regions of (R, r), with taking only the wave 
function in the breakup channels. The expression of Eq. (|14|) is unphysical, because 
(i) breakup channels can contribute not only the incomplete but also complete fusion 
processes, (ii) a possible contribution from the elastic channel is naively disregarded, 
and (hi) couplings between the channels of the three-body system, which have been 
included in the calculation of \P(R, r), are neglected in the evaluation of Ojp v ; the 
expression of the complete fusion given in Ref. 4), o"cF,prev, has similar issues. 

We apply the new method for calculating complete and incomplete fusion cross 
sections to the deuteron induced reactions on the 7 Li target for 10 MeV < < 
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50 MeV, where E^ is the deuteron incident energy in the laboratory system. We 
use the CDCC codes cdcdeu and hicadeu 15 ) to obtain the CDCC wave function 
with assuming intrinsic spins of p, n, and Li to be zero. As for V pn , we adopt 
the Ohmura potential 16 ) that reproduces the deuteron energy in the ground state, 
i.e., Eq = —2.23 MeV. In the calculation of the p-n discretized continuum states, 
we include the s- and d-wave states; the maximum relative wave number /c max is 
determined by 

U = \^2u r {E™ - M) 

with E9 M being the d- 7 Li relative energy. The p-n continuum state is divided into 
4 bin states, for each of the s- and d-waves. As for the p- 7 Li and n- 7 Li optical 
potentials, we use the parameter sets in Ref. 17) except that the spin-orbit terms 
are neglected in this study. 

In the present formalism, the absorption radius r'^° is assumed to be a free 
parameter. In fact, it is found that the results of the incomplete fusion cross sections 
calculated with rjf 3 = 4 and 5 fm differ from each other by about 30%. Therefore, 
in this study, we determine r^ b at E\ = 40 MeV so that the a[p agrees with the 
result of the Glauber model calculation; the latter, added by the elastic breakup 
contribution calculated with CDCC, reproduces the experimental DDX data very 
well at the same incident energy. The absorption radius thus determined is 4.0 fm, 
which is used for both p and n in all calculations in this study. 

Figure [3] shows the results calculated with the new method; the dash-double- 

(p) (n) 

dotted, solid, dashed, and dash-dotted lines represent cttf, ccf> <Tjf , and <r IF , 
respectively. One sees that <tcf has the largest contribution to cttf in the energy 
region of our interest. Another important feature is that the energy dependence of 
a^} is significantly different from that of ; at E\ = 10 MeV, a[p is three times as 

(n) 

large as CTj F . This is due to the difference in the energy dependence of W p and W n , 

1000 
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b 400 
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10 20 30 40 50 
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Fig. 3. (Color online) Results of cttf (dash-double-dotted line), ctcf (solid line), a jp (dashed line), 
and (Tjp (dash-dotted line) for the deuteron induced reaction on 7 Li as a function of the incident 
energy E^. The dotted line represents the elastic breakup cross sections calculated with CDCC. 
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method, compared with cr^p v (squares) and a^' prev (triangles) following the previous defini- 
tion given in Refs. 4) and 13). b) Comparison between the results of the complete fusion cross 
sections calculated with the present (solid line) and previous (dots) method. 



Fig. 4. (Color online) a) a^ 1 (dashed line) and <7jp' (dash-dotted line) calculated with the present 



» 



i.e., W p (W„) at low energy is more (less) absorptive than that at around 40 MeV. 
The elastic breakup cross section <teb is also shown by the dotted-line in Fig. [3l 
The total neutron emission cross section, except for those by the compound and 

preequilibrium processes, can be evaluated as the sum of crj^ and cteb- It is found 

(p) 

that contribution of <teb is much smaller than that of <7j F , which is consistent with 
the conclusion of Ref. 11). On the other hand, if we consider the proton emission 
cross section below 20 MeV, which consists of a^) and <teb> the two contributions 
are comparable. 

Next we show in Fig. [J] the incomplete and complete fusion cross sections calcu- 
lated with the previous expressions in Refs. 4) and 13), compared with the results 

of the present study. The squares and triangles in the left panel show, respectively 

(p) (n) 

°IF prcv an d °iF P rev> an( ^ *he d°ts in ^ e ri ght panel show 0"cF,prev ^he li nes shown 
in the panels are the same as in Fig. [3l One sees clearly that the previous prescrip- 
tion gives much smaller (larger) incomplete (complete) cross sections than those 

(p) 

obtained by the present calculation. Our result of afp at 40 MeV, by definition, can 
be interpreted as an experimental value. Thus, the prescription given in Refs. 4) 
and 13) do not work at all at least for the inclusive 7 Li(<i, n) reaction concerned. In 
other words, Fig. H] clearly shows the importance of including elastic channel in the 
evaluation of crjp (c = p or n) as in Eqs. (|12p and (|13p . Similarly, inclusion of the 
breakup channels in the calculation of acF is also important. 

As mentioned above, in Ref. 11), the contribution of the proton stripping process 
to the inclusive (d, n) cross section at 40 MeV was shown to be described very well 
by the Glauber model. It is thus interesting whether the Glauber model calculation 
works or not at even lower energies. For this purpose, in the left-upper panel in 
Fig. [5l we compare o"j F (short-dashed line) calculated with the present method, 
with OgTR (squares) obtained by the Glauber model calculation. 11 ) Surprisingly, 
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Fig. 5. (Color online) a) Comparison of Ujp (short-dashed line) and tjjp' (dash-dotted line) ob- 
tained by the present calculation with <t S tr (squares) and o" S x R (triangles) by the Glauber 
model. 11 ' b) Complete fusion cross sections calculated with the Glauber model (dots) and the 
present method (solid line), c) The elastic breakup cross section calculated with CDCC (dot- 
ted line) compared with that with the Glauber model (inverse triangles), d) The total reaction 
cross section obtained by CDCC (long-dashed line) is compared with that by the Glauber model 
(diamonds). 



the two results agree well each other not only above 40 MeV but also at low energies 
down to 10 MeV. This is also the case with the neutron stripping process; one sees 
the good agreement between a^p (dash-dotted line) and croTR (triangles). Note that 
we use the absorption radius of 4.0 fm determined at 40 MeV in all the calculation 
as mentioned above. Thus, we conclude that the Glauber model calculation for the 
incomplete fusion cross sections is expected to work even at lower energies down to 10 
MeV. On the other hand, as shown in the left-lower panel, the results of the complete 
fusion cross section ctq F obtained by the Glauber model (dots) 14 ) significantly deviate 
from those obtained by the present study, i.e., <tcf (solid line); even the energy 
dependence is different. 

A possible reason for the success of the Glauber model in describing stripping 
processes is, as discussed in Ref. 11), that the contribution from the nuclear surface 
region is dominant, where the depth of the optical potential is so shallow that the 
Glauber model works well. On the other hand, since the complete fusion process 
takes place in the nuclear interior, the Glauber model does not work even at 50 
MeV. It is numerically confirmed that when we make artificially the nucleon- 7 Li 
optical potential shallow, the difference shown in the left-lower panel becomes small, 
while the features of the incomplete fusion cross sections (left-upper panel) have no 
changes. 

In the right-upper panel, we show the results of the elastic breakup cross section 
calculated with CDCC (dotted line) and the Glauber model (inverse triangles). The 
difference between the two is quite small above 10 MeV, which is found to be mainly 
due to the adiabatic approximation used in the Glauber model. The quite big dif- 
ference at 10 MeV will come from the invalidity of the eikonal approximation also 
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assumed in the Glauber model. The total reaction cross section, which is the sum 
of the total fusion cross section and the elastic breakup one, calculated with CDCC 
(the Glauber model) is shown by the long-dashed line (diamonds) in the right-lower 
panel. The main part of the difference of the two comes from that in the complete 
fusion cross sections. 

In summary, we propose a new method for evaluating the complete and incom- 
plete fusion cross sections separately by means of CDCC. The separation of the two 
is carried out by the physics condition on the absorption for each fusion process. The 
absorption radius included in the present formalism is determined using the result of 
the proton stripping cross section for the 7 Li(d, n) reaction at 40 MeV calculated by 
the Glauber model, which was shown to reproduce the corresponding experimental 
value. The new method is applied to the 7 Li(d, n) reaction from 10 to 50 MeV. The 
complete fusion cross section is found to have the largest contribution to the total 
fusion cross section. The (d,p) and (d,n) incomplete fusion cross sections show quite 
different energy dependence, because of that in the imaginary parts of the p- 7 Li and 
n- Li optical potentials. It is found that in the all energy region considered, the (d,p) 
and (d, n) incomplete fusion cross sections obtained by the Glauber model agree well 
with those obtained by the present calculation with CDCC. On the other hand, 
the two model calculations give significantly different results of the complete fusion 
cross section, even at 50 MeV. The complete and incomplete fusion cross sections 
obtained by the previous method of Refs. 4) and 13) are found to be inaccurate. 
Extension of the present framework to calculate the DDX is a very important future 
work. A method to divide the complete and incomplete processes unambiguously, 
i.e., without the absorption radius, will also be desirable. 

We would like to thank Y. Watanabe and T. Ye for fruitful discussions and pro- 
viding the numerical results. We also acknowledge helpful discussions with K. Hagino. 
SH thanks Y. Aoki for valuable discussions. 
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